Giant dipole resonance with exact treatment of thermal fluctuations 
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The shape fluctuations due to thermal effects in the giant dipole resonance (GDR) observables 
are calculated using the exact free energies evaluated at fixed spin and temperature. The results 
obtained are compared with Landau theory calculations done by parameterizing the free energy. 
The Landau theory is found to be insufficient when the shell effects are dominating. 
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The study of structural transitions as a function of both angular momentum and temperature has been one among 
the fascinating aspects of highly excited nuclei in recent years 0, IS 3 ■ The Giant Dipole Resonance (GDR) 
' studies have been proved to be a powerful tool to study such hot and rotating nuclei|al and recently the domain of 
fSJ GDR spreads rapidly over different areas of theoretical and experimental interest 0, UH 1^ • The GDR observations 
(— ( I provide us information about the geometry as well as the dynamics of nuclei even at extreme limits of temperature 
I (T), spin (J) and isospin (r)- In the past most of the GDR measurements in hot nuclei were made at moderate and 
, high T |H S EL 01 S liol. Ej • Several experiments have been carried out recently to study the GDR states at low 
^sj ■ temperatures [l|, 01 13 ■ Hence the theories which were successful in the high T regime should be now scrutinized with 
the low temperature observations as well as with theories incorporating properly the microscopic effects (such as shell 
effects) which are dominant at low temperatures. While dealing with the thermal shape fluctuations, free energy 
parameterizations such as Landau theory [T^ll3llT^ are usually employed to do timesaving calculations. In this work 
we survey the applicability of Landau theory by demanding consistency with exact calculations done without any 
' parameter fitting. 

The theoretical approach we follow is of three fold with models for 1) shape calculations, 2) relating the shapes to 
T— I GDR observables and 3) considering the shape fluctuations due to thermal effects. For shape calculations we follow 
the Nilsson-Strutinsky (NS) method extended to high spin and temperature p^lTslIT^ . The total free energy (Ftot) 
^•^^ at fixed deformation is calculated using the expression 

; -FtOT = -BlDM + ^ <5F + itj(/ciassical ■ (1) 

o ' . . . |— I . 

^ ' The liquid-drop energy (£'ldm) is calculated by summing up the Coulomb and surface energies 115^ corresponding to 
(~| , a triaxially deformed shape defined by the deformation parameters (3 and 7. The classical part of spin (/classical) is 
^ ' obtained from the rigid-body moment of inertia with surface diffuseness correction |15| . The shell correction {5F) is 
. J—^ , the difference between the deformation energies evaluated with a discrete single-particle spectrum and by smoothing 

(averaging) that spectrum {5F — F—F). Similarly the shell correction corresponding to the spin is given hy 51 = /— /. 
^ , To calculate the shell corrections for energy and spin, we use the triaxially deformed Nilsson model together with the 
- - I Strutinsky's prescription. The single-particle energies (e^) and spin projections {rrii) are obtained by diagonalizing 

the triaxial Nilsson Hamiltonian in cylindrical representation upto first twelve major shells. At finite temperatures 

the free energy is given by 

F = ^e,n, -T^s,, (2) 

i=l 1=1 

where Si are the single-particle entropy and rii are the occupation numbers which follow Fermi-Dirac distribution 
given by 

ni — ■ (3) 

H-exp(2Y^) ' 
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The chemical potential A is obtained using the constraint X^i^i '^i ~ where N is the total number of particles. 
The total entropy S — X^i^i '^^'^ represented in terms of occupation numbers as 

5 = -^Kln?i, -(l-7i,)ln(l-n0] ■ (4) 

i=l 

The shell correction energy at finite temperature is calculated using the expression 0, 0| 

F = ^e^fii ~ T^Si , (5) 

i i 

where ni = f{^) n-iix) x , and Si — fi^) Si{x) x. The integrals appearing in Eq.jSJ) are evaluated numerically 
using the Hermite-Gauss quadrature. Apart from the numerical evaluation of integrals, this method gives exact 
temperature dependent shell corrections. For the spin distribution, the Strutinsky smoothed spin can be derived in a 
similar way leading to the expression / = "^i^j ^-iid hence the shell correction for spin is 

oo oc 

5/ = ^ mim - ^ mjHi . (6) 

1=1 i=l 

The most probable shapes are obtained by minimizing the total free energy (Q. 

The nuclear shapes are related to the GDR observables using a model Iq comprising an anisotropic harmonic 
oscillator potential with separable dipole-dipole interaction. In this formalism the GDR frequencies in laboratory 
frame are obtained as 

= (l + r/)i/V , (7) 
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(8) 



+8n'{i + v){Lol + ujl)]^Y Tn, (9) 

where n is the cranking frequency, lOx, ujy, i^z are the oscillator frequencies derived from the deformation of the 
nucleus and 77 is a parameter that characterizes the isovector component of the neutron and proton average field. The 
GDR cross sections are constructed as a sum of Lorentzians given by 

= E , (10) 



where Lorentz parameters -Em, Cm and F are the resonance energy, peak cross-section and full width at half maximum 
respectively. Here i represents the number of components of the GDR and is determined from the shape of the nucleus 
[isL [13 • It is to be noted that these Lorentz lines are non- interfering, but F^ is assumed to depend on energy. The 
energy dependence of the GDR width can be approximated by 

F, « OmGEl^ . (11) 

The peak cross section am is given by 

2 NZ 1 

a,n = 60 ^-0.86(1 + a). (12) 



The parameter a which takes care of the sum rule is fixed at 0.3 for all the nuclei considered in this work. This 
parameter has more effect on the peak cross section. In most of the cases we normalize the peak with the experimental 
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data and hence the choice of a has less effect on the results. The other parameter r] varies with nucleus so that the 
ground state GDR centroid energy is reproduced. The choice for ^°Zr, ^^Mo is ij = 2.6 and for ^"^Pb it is 77 = 3.4. 
For calculating the GDR width, only the power law 1)11(1 is used in this work and no ground state width is assumed. 

The relation between nuclear shape and GDR cross section is not straightforward especially in hot nuclei where 
large-amplitude thermal fluctuations of the nuclear shape play an important role When the nucleus is observed 
at finite excitation energy, the effective GDR cross-sections carry information on the relative time scales for shape 
rearrangements. Hence in the case of hot nuclei, for a meaningful comparison of experimental and theoretical values, 
the thermal shape fluctuations should be taken care properly. In the case of hot and rotating nuclei there can be 
fluctuations in the orientation of the nuclear symmetry axis with respect to the rotation axis. The general expression 
for the expectation value of an observable O incorporating both thermal and orientation fluctuations is given by 

/P[a]e-^(™.")/^(^.T.^)-3/20 

{U)l3,j,n J p[c,]e-F(TJ;/3,7,n)/T(^. 1.^)^3/2 ' ^^'') 

where il = (0, 9, "0) are the Euler angles specifying the intrinsic orientation of the system, uj-X-uj = I^'x' cos^ cj) sin^ 9 + 
ly'y' sin^ (j) sin^ 9-\-Iz'z' cos^ 6 is the moment of inertia about the rotation axis w given in terms of the principal moments 
of inertia Ix'x' : ly'y' 1 ly'x', and the volume element 'D[a.] = sin37| d(3 sm9 d9 dcj). 

The study of thermal fluctuations by numerical evaluation of Eq. I|13fl in general requires an exploration of five 
dimensional space spanned by the deformation and orientation degrees of freedom, in which a large number of 
points are required in order to assure sufficient accuracy (especially at finite angular momentum). Hence certain 
parameterizations were developed [2 ^3 represent the free energy using functions that mimic the behaviour of the 
NS calculation as closely as possible. One such parameterization is the Landau theory of phase transitions, developed 
by Alhassid et al [l^ . Here the free energy is expanded in terms of certain temperature dependent constants which 
are to be extracted by fitting with the free energy calculations at fixed temperatures from the NS method. Moreover, 
once the fits involving free energy and moment of inertia are made for the non-rotating case, the calculations can be 
extended to higher spins using the relation jl^ 

F(T, /; 7, n) = F{T, = 0; 7) + ■ (14) 

Hence this theory offers an economic parameterization to study the hot rotating nuclei. We have employed Landau 
theory in its extended form as given in Refs. dill. 

With recent computing facilities, it is possible to perform the thermal fluctuation calculations exactly by computing 
the integrations in Eq. (|13|l numerically with the free energies and the observables being calculated exactly at the 
integration (mesh) points. In this way the calculations can be done more accurately without using any parameteri- 
zation and consequent fltting. In this work we have performed such calculations, however, neglecting the orientation 
fluctuations. This enables us to perform the integration in the deformation space only which at present is two dimen- 
sional having the deformation parameters (3 and 7. A 32 point Gaussian quadrature has been used to evaluate the 
integrals numerically. While performing fluctuation calculations in this way, the free energy, GDR cross section and 
width at any given spin is obtained by tuning the cranking frequency to get the desired spin. It has to be noted that 
the NS free energy calculations at mesh points (32 x 32) have to be done only once for a fixed spin and temperature 
while calculating the GDR cross sections at fixed spin and temperature. 

Now we compare our calculated results obtained by using the extended Landau theory and the exact method. In 
Fig. 1 as the representative cases we show for the hot rotating ^°Zr and ^^Mo nuclei, the results of our GDR cross 
section calculations along with the experimental results d|. At experimentally observed temperatures, the most 
probable shapes are found to be /3=0.0, 0.09 and 0.1 for the spins 9ft, 22h and 33ft respectively, with 7 = —180° for 
all the spins. The thermally averaged value of the deformation, /3 is found to be 0.29, 0.33 and 0.41 for the spins 9ft, 
22ft and 33ft respectively, and the averaged 7 is found to be 7 w 148°. These values correspond to Landau theory 
calculations and are consistent with those reported earlier . As similar to previous observations (22| it is clear from 
Fig. 1 that thermal fluctuations play crucial role whereas the orientation fluctuations are negligible while calculating 
the GDR cross sections. Also this justifies the omission of orientation fluctuations in our exact method. Importantly 
we find that the results of the thermal fluctuation calculations using the macroscopic extended Landau theory and 
the present exact approach, are exactly similar (they exactly overlap in left panel of Fig. 1) at the experimentally 
observed temperatures. This equivalence is in spite of the differences in the free energies, from NS calculations and 
Landau theory, at larger deformations. The thermal fluctuations are strong enough to make the GDR cross sections 
insensitive to minor changes in the free energy surfaces. We examined the situation at lower temperatures and the 
cross sections calculated using two methods at T = 0.5 MeV are shown in the right panel of Fig. 1. To enable better 
comparison we show results without orientation fluctuations. It is evident that the Landau theory results deviate at 
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FIG. 1: GDR cross sections for the nuclei ^"Zt and ^^Mo. Left: Experimental data represented by solid squares are taken from 
ref. The solid lines represent calculations with orientation and thermal fluctuations, dashed lines correspond to thermal 

fluctuations alone and the dash-dotted line correspond to most probable shapes. Right: Results obtained at T = 0.5 MeV with 
Landau theory (dashed lines) and exact method (solid line). 

/ = 33ft. This can be ascribed to the spin driven shell effects. The rapid change in the single-particle level structure 
with spin is well known and this may lead to sharp changes in the shell corrections. The Landau theory could not 
account for this as the free energies at higher spins are obtained from the free energies and moment of inertia at / = 
using Eq. (fTIjl . 

In Fig. 2 we present the calculated GDR widths of ^°*Pb at / = Oft along with experimental results. The calculations 
are performed with 1) the liquid-drop model (LDM) free energies and Landau theory, 2) NS free energies and Landau 
theory and 3) NS free energies with exact treatment of fluctuations. The strong shell corrections for spherical shape 
results in large difference in the deformation energies between spherical and deformed configurations. This leads 
to attenuation of thermal fluctuations at lower temperature and hence the obtained widths are much lower when 
compared to liquid-drop model results. The magnitude of this attenuation comes out to be different in methods 2 
and 3. In the presence of strong shell effects at lower temperatures, the Landau theory results deviate considerably 
from the results of exact calculations. Above T ~ l.SMeV the results from two methods are the same. The deviations 
from the NS calculations in the free energy obtained by Landau theory at larger deformations are not reflected 
in the GDR observables calculated at higher temperatures due to the thermal fluctuations and the weakening of shell 
effects. The results of our calculations suggest that the parameterization of Landau theory is good enough to explain 
the GDR properties of hot rotating nuclei in the absence of strong shell effects. The shell effects can be treated in a 
better way using the parameterization suggested in Ref. However, spin driven shell effects cannot be explained 

in such a formalism also as the difference in energy due to rotation comes just through the moment of inertia. 

To summarize, in this work the thermal fluctuations are dealt in an exact way without any parameter fitting. We 
have carried out a case study of the GDR properties in the nuclei ^°Zr, ^^Mo and ^'^^Pb and our results are well 
in conformity with experimental results. Comparison of our present approach with the thermal fluctuation model 
comprising Landau theory has been brought out. The Landau theory is found to be insufficient to explain GDR 
properties in the presence of strong shell effects. 
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FIG. 2: GDR width in ^"^^Pb. The results obtained using liquid drop model (dash-dotted line), Landau theory (dashed line 
and the exact calculations (solid line) are compared. Experimental data represented by solid squares are taken from ref. [ij 
and the revised data j.24j are represented by open circles. 
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